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Abstract: We investigate the combined influence of chemical reaction and radiation absorption on mixed convective flow in 
a vertical channel with oscillatory wall temperature and concentration . The non-linear coupled partial differential equations 
governing the flow heat and mass transfer are solved by a perturbation technique. The effect of various forces acting on the 
fluid system is analyzed by graphical representation of the velocity, temperature and concentration. 
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I. INTRODUCTION 

Combined heat and mass transfer problems with chemical reaction are of importance in many processes and have, 
therefore, received a considerable amount attention in recent years. In processes such as drying, evaporation at the surface of 
a water body, energy transfer in a wet cooling tower and the flow in a desert cooler, heat and mass transfer occur 
simultaneously. 

We are particularly interested in cases in which diffusion and chemical reaction occur at roughly the same speed. 
When diffusion is much faster than chemical reaction, then only chemical factors influence the chemical reaction rate; when 
diffusion is not much faster than reaction, the diffusion and kinetics interact to produce very different effects. The study of 
heat generation or absorption effects in moving fluids is important in view of several physical problems, such as fluids 
undergoing exothermic or endothermic chemical reaction. Due to the fast growth of electronic technology, effective cooling 
of electronic equipment has become warranted and cooling of electronic equipment ranges from individual transistors to 
main frame computers and from energy suppliers to telephone switch boards and thermal diffusion effect has been utilized 
for isotopes separation in the mixture between gases with very light molecular weight (hydrogen and helium ) and medium 
molecular weight. 

Muthucumaraswamy and Ganesan (12) studied effect of the chemical reaction and injection on flow characteristics 
in an unsteady upward motion of an unsteady upward motion of an isothermal plate. Deka et at. (4) studied the effect of the 
first order homogeneous chemical reaction on the process of an unsteady flow past an infinite vertical plate with a constant 
beat and mass transfer. Chamkha (3) studies the MHD flow of a numerical of uniformly stretched vertical permeable surface 
in the presence of heat generation/absorption and a chemical reaction. The effect of foreign mass on the free-convection flow 
past a semi-infinite vertical plate were studied (5) Chamkha (3) assumed that the plate is embedded in a uniform porous 
medium and moves with a constant velocity in the flow direction in the presence of a transverse magnetic field. Raptis and 
Perdikis (19) studied the unsteady free convection flow of water near 4 C in the laminar boundary layer over a vertical 
moving porous plate. 

In spite of all these studies, the unsteady MHD free convection heat and mass transfer for a heat generating fluid 
with radiation absorption has received little attention. Hence, the main objective of the present investigation is to study the 
effects of radiation absorption, mass diffusion, chemical reaction and hear source parameter of hear generating fluid. The 
process of free convection as a mode of heat transfer has wide applications in the fields of Chemical Engineering, 
Aeronautical and Nuclear power generation. It was shown by Gill and Casal (6) that the buoyancy significantly affects the 
flow of low Prandtl number fluids which is highly sensitive to gravitational force and the extent to which the buoyancy force 
influences a forced flow is a topic of interest. Free convection flows between two long vertical plates have been studied for 
many years because of their engineering applications in the fields of nuclear reactors, heat exchangers, cooling appliances in 
electronic instruments. These flows were studied by assuming the plates at two different constant temperatures or 
temperature of the plates varying linearly along the plates etc. The study of fully developed free convection flow between 
two parallel plates at constant temperature was initiated by Ostrach (16). Combined natural and forced convection laminar 
flow with linear wall temperature profile was also studied by Ostrach (17). The first exact solution for free convection in a 
vertical parallel plate channel with asymmetric heating for a fluid of constant properties was presented by Anug (1). Many of 
the early works on free convection flows in open channels have been reviewed by Manca et al. (11). Recently, Campo et al. 

(2) considered natural convection for heated iso-flux boundaries of the channel containing a low-Prandtl number 
fluid. Pantokratoras (18) studied the fully developed free convection flow between two asymmetrically heated vertical 
parallel plates for a fluid of varying thermophysical properties. However, all the above studies are restricted to fully 
developed steady state flows. Very few papers deal with unsteady flow situations in vertical parallel plate channels. 

Transient free convection flow between two long vertical parallel plates maintained at constant but unequal 
temperatures was studied by Singh et al.(20). Jha et al. (9) extended the problem to consider symmetric heating of the 
channel walls. Narahari et al. (15) analyzed the transient free convection flow between two long vertical parallel plates with 
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constant heat flux at one boundary, the other being maintained at a constant temperature. Singh and Paul (20) presented an 
analysis of the transient free convective flow of a viscous incompressible fluid between two parallel vertical walls occurring 
as a result of asymmetric heating / cooling of the walls. Narahari (14) presented an exact solution to the problem of unsteady 
free convective flow of a viscous incompressible fluid between two long vertical parallel plates with the plate temperature 
linearly varying with time at one boundary, the other boundary being held at constant. There are many reasons for the flow to 
become unsteady. When the current is periodic due to on-off control mechanisms or due to partially rectified a-c voltage, 
there exist periodic heat inputs. Hence, it is important to study the effects of periodic heat flux on the unsteady natural 
convection, imposed on one of the plates of a channel formed by two long vertical parallel plates, the other being held at a 
constant initial fluid temperature. Recently Narahari(15) has discussed the unsteady free convection flow of dissipative 
viscous incompressible fluid between two long vertical parallel plates in which the temperature of one of the plates is 
oscillatory whereas that of the other plate is uniform. Haritha (7) has analysed unsteady convective heat transfer of 
dissipative viscous fluid through a porous medium confined in a vertical channel on whose walls an oscillatory temperature 
is prescribed. Ibrahim et al. (8) have studied the effect of chemical reaction and radiation absorption on the unsteady MHD 
free convection flow past a semi infinite vertical permeable moving plate with heat source and suction Kesavaiah et al (10) 
have studied the effect of the chemical reaction and radiation absorption on an unsteady MHD convective Heat and Mass 
Transfer flow past a semi-infinite vertical permeable moving plate embedded in a porous medium with heat source and 
suction 

In this paper we analyse the effect of chemical reaction and radiation absorption on unsteady convective Heat and 
Mass Transfer flow of a viscous fluid in a Vertical channel on whose walls an oscillatory temperature is prescribed. 
Approximate solutions to coupled non-linear partial differential equations governing the flow, heat and mass transfer are 
solved by a perturbation technique. The velocity, temperature, skin friction, concentration and rate of heat and mass transfer 
are discussed for different variations of Sc, N, K, Q[ y. 




II. FORMULATION OF THE PROBLEM 

We consider the flow of a viscous incompressible chemically reacting fluid in a vertical channel bounded by flat 
walls in the presence of constant heat sources. We choose a Cartesian coordinate system 0(x y) with walls at y = ± 1 by using 
Boussinesq approximation we consider the density variation only on the buoyancy term. The equation governing the flow to 
heat and mass transfer are 



Equation of Linear Momentum 
du U d 2 U _ 
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Equation of State 

P- Po =~J3 0 (T~T 0 )-J3*(C- C 0 ) (4) 

where u is the velocity component in x-direction, T is the temperature, p is the pressure, p is the density, a i 



(he 



electrically conductivity, |i e is the magnetic permeability, k is the coefficient of porous permeability, |i is dynamic viscocity, 
k f is coefficient of thermal conductivity p is the coefficient of volume expansion, Q is the strength of heat source, P* is the 
volumetric coefficient of expansion with mass fraction, D i is the chemical molecular diffusivity and K' is the coefficient of 
chemical reaction. Qj is radiation absorption parameter. 
The boundary conditions are 
u = 0,T = Ti ,C=Clat y = -L 
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= 0, T = T, + e(T 2 - Ti) cos( rot), C = d + e(C 2 - Cj) cos( rot)at y = +L (5) 



On introducing the non-dimensional variables. 

u _ T-T, 

u = , y' = y/L, 9 = L , 

vIL T 7 -T t 



_ (c-c t ) 
~(c 2 - Cl ) 



Equ; 



>ns (2.1) - (2.3) reduce to (dropping the dashes) 
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(Grashoff number) 
(Prandtl number) 
(Heat source parameter) 
(Schmidt number) 
(Wormsly number) 
(Buoyancy ratio) 
(Chemical reaction parameter) 
(radiation absorption parameter) 



The transformed boundary conditions are 
u = 0, 9 = 0, <j)= 0 at y = -1 
u =0, 9 = 1 + e cos(rot), c = 1 + e cos(rot) at y = +1 



(9) 



III. METHOD OF SOLUTION 

In view of the boundary conditions (5) we assume 
u = u 0 + e e" U| -\ 

9 = 9 0 +ee"9i (10) 
<t> = <t>o +e e"<|)i L 

substituting (10) in the equatinns6-8 and comparing harmonic & Non harmonic terms we get 

d 2 u Q 

w 

^4"- CD 1 +//")//, = -G(0, +Nfi) (12) 

dy 
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d% 
By 2 ' 
d\ 
dy 2 
d\ 
dy 2 



-iPy 2 0 l + Q l </> l = 0 (14) 
-K^ 0 =0 (15) 
-(K +iScy 2 )0 1 = 0 (16) 



III. SOLUTIONS OF THE PROBLEM 

The solutions of (1 1)- (16) subject to the boundary conditions are 
_ Chfiy | Shfty 
IChp, 2Shp, 
8 0 = b, y 2 + b 2 Ch P,y + b 3 Sh pi y + a 3 y + a 4 

u 0 = b 4 y 4 + b 7 y 3 + b 8 y 2 + a 5 y + a 6 + (b 5 + b 9 ) Ch Pi y) + (b 6 + b 10 ) Sh Pjy. 

_ Chfry Shfry 
1 2Chj3 3 2Shj3 3 

9, = k 49 Chp 3 y + k 50 Sh p 3 y + k 5 , Ch p 2 y + k 52 Sh p 2 y 

ui = k 55 Ch p 3 y + k 56 Sh p 3 y + k 57 Chp 2 y + k 58 Shp 2 y + k 59 Ch p 4 y + k 60 Sh p 4 y. 
Where J3 2 = a fi 2 = K fi 2 = iy 2 P J3 2 = iy 2 Sc 

IV. NUSSELT NUMBER and SHERWOOD NUMBER 

The rate of heat transfer (Nusselt number) at the walls y = ±1 is given by 



and the corresponding expressions are 

(Nu) y= _, = a 21 + Ec[a 2S + a 29 + a^Shft + a 3] Chp 1 ] + (0. 01). E 33 .a 32 

(Nu) Y=+l = a 2l + Ec[a 22 + a 23 + a 2A Shp, + a 25 Chp l \ + (0. 01). E 33 .a 26 
The rate of mass transfer (Sherwood number) at the walls y = ±1 is given by 

And the corresponding expressions are 

(S/z) v= _, =a 35 + (0.01)£33.a 36 

(Sh) y=+ , =a 33 + (0.0l)E 33 .a 34 
where a b a 2 , a 36 are constants shown in appendix. 

V. RESULTS AND DISCUSSION 

In this analysis we investigate the effect of chemical reaction on mixed convection heat and mass transfer flow of a 
viscous fluid through a porous medium in a vertical channel on whose walls oscillatory temperature and concentration are 
prescribed. 

The velocity (u) is shown in figures 1-4 for different values of N,Sc , K, The variation of u with buoyancy ratio 
N shows that when the molecular buoyancy force dominates over the thermal buoyancy force, the velocity enhances in the 
left half and reduces in the right half when the buoyancy forces act in the same direction and for the forces acting in 
opposite directions u reduces in the left half and enhances in the right half(fig.l). Fig 2 represents the variation of 'u' with 
Sc. Lesser the molecular diffusivity smaller 'u' in the flow region, and for further lowering of the diffusivity the velocity 
enhances in the left half and reduces in the right half and for still lowering of the molecular diffusivity the velocity enhances 
in the entire flow region. The effect of chemical reaction on u is shown in fig 3. It is found that an increase in k < 1.5 
enhances lul in entire region and for higher k > 2.5, lul enhances in the left half and reduces in the right half. The variation of 
'u' with radiation absorption parameter 'Ch' is shown in fig 4. Fixing the other parameters. It is found that an increase in Qj 
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< 2 leads to an enhancement in lul and for further higher Q, = 4, lul reduces in the flow region and for still higher Qi > 6, lul 
reduces in the left half and enhances in the right half. 




The Non-dimensional temperature distribution (9) is shown in fig 5- 8 for different parametric values N,Sc , K, Q[ . 
The variation of 9 with the buoyancy ratio N shows that the actual temperature experiences a depreciation with increasing in 
|N| irrespective of the directions of the buoyancy forces(fig.5). The variation of '9' with Schmidt number Sc shows that 
lesser the molecular diffusivity larger the actual temperature in the left half and smaller in the right half, and for further 
lowering of the molecular diffusivity smaller the actual temperature in the entire flow region and for still lowering of the 
molecular diffusivity smaller the actual temperature in the left half and larger in the right half (fig 6). From fig 7 we notice 
that for smaller and larger values of the chemical reaction parameter K, the actual temperature depreciates in the flow region 
and for any intermediate value of K = 2.5, the actual temperature experiences an enhancement in the flow region. The 
variation of '9' with radiation absorption; parameter Qi shows that the actual temperature depreciates appreciabelly with 
increase in Qj in the entire flow region (fig.8). 
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The non-dimensional concentration 'C is shown in fig 9-11 for different values of Sc, K and y. From fig 9 we 
notice that lesser the molecular diffusivity larger the actual concentration in the flow region. An increase in the chemical 
reaction parameter 'K' results in a depreciation in the concentration in the entire flow region (fig. 10). An increase in the 
Wormsely number (y) leads to an enhancement in the actual concentration in the flow field(fig.l 1). 



i ii m iv 

Sc Cl24 as L3 XOl 



The rate of heat transfer at y = ±1 are exhibited in tables 1-4 for different parametric values. It is found that the rate 
of heat transfer increases with increase in IG| (or) M (or) K, thus higher the Lorentz force larger the Nusselt number at both 
the walls. An increase in the radiation absorption parameter Qi enhances INul at y=+l and reduces at y = -1 (tables 1 and 3). 

From tables 2 and 4 we find that the rate of heat transfer depreciates at y = +1 and enhances at y = -1. The variation 
of Nu with heat source parameter a shows that the rate of heat transfer reduces at y = +1 and enhances at y = -1 with 
increase in a > 0, while an increase in a < 0 reduces INul at y = +1. 

The rate of mass transfer at y = ±1 is shown tables 5-6 for different Sc, K and y. The rate of mass transfer enhances 
with Schmidt number Sc. Thus higher the molecular diffusivity larger IShl at y = ±1. The variation of Sh with chemical 
reaction parameter K shows that the rate of mass transfer enhances with K at y = +1 while at y = -1 it enhances with K < 1.5 
and reduces with higher K > 2.5. With reference to variation of Sh with Wormsely number (y) exhibits that IShl reduces with 
increase in y < 4 and enhances with y > 6, while at y = -1 it experiences an enhancement with y(tables.5-7). 
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Table 1 

Nussclt Number (Nu) at y = + 1 



10.62459 53.76271 



2.29075 2.19075 



21.04668 
2.15887 
21.25105 



97.73100 485.93610 40.50675 
10.65264 53.77809 4.45838 
97.81516 485.9823 40.63005 



76.52689 
8.63823 



1124.36 
125.4219 
1 125.023 



2.25887 2.20887 
22.25105 21.95105 



20.00887 
20A5105 



It Number (Nu) at y = - 1 



24.12496 38.2541 
35.0605 43.4371 



19.077 



45.25756 54.98854 
66.0725 77.4422 



65.200 
86.156 



2.48786 
6.46406 



1.25672 
1.63609 



Nusselt Number (Nu) at y = 



Sifmroud Numbf r (Sh) at 
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Table 6 



0.24 
0.6 



1.44945 
1.71269 



1.52393 
1.68234 



1.05223 
1.16664 



0.85282 
0.93067 



1.50818 
1.90167 



2.49153 
3.24322 



-101.3 
2.01 



1.98342 
2.15846 



1.88206 
2.02504 



1.32729 
1.44943 



1.05200 
1.15108 



2.42227 
" 2.76520 



1.03263 
1.55103 



VL CONCLUSION 

An attempt has been made to investigate the combined influence of chemical reaction and radiation absoiption on 
the unsteady convective heat and mass transfer flow in a vertical channel technique using a regular perturbation technique 
the non-linear coupled equations has been solved. The important conclusions of this analysis are 

1. Lesser the molecular diffusivity smaller 'u' in the flow region, and for further lowering of the diffusivity the velocity 
enhances in the left half and reduces in the right half and for still lowering of the molecular diffusivity the velocity 
enhances in the entire flow region. An increase in k < 1.5 enhances lul in entire region and for higher k > 2.5, lul 
enhances in the left half and reduces in the right half.An increase in Qi < 2 leads to an enhancement in lul and for further 
higher Qi = 4, lul reduces in the flow region and for still higher Q { > 6, lul reduces in the left half and enhances in the 
right half. 

2. The actual temperature experiences a depreciation with increasing in |N| irrespective of the directions of the buoyancy 
forces. Lesser the molecular diffusivity larger the actual temperature in the left half and smaller in the right half, and for 
further lowering of the molecular diffusivity smaller the actual temperature in the entire flow region and for still 
lowering of the molecular diffusivity smaller the actual temperature in the left half and larger in the right half. 

3. For smaller and larger values of the chemical reaction parameter K, the actual temperature depreciates in the flow region 
and for any intermediate value of K = 2.5, the actual temperature experiences an enhancement in the flow region. The 
actual temperature depreciates appreciabelly with increase in Q] in the entire flow region. 

4. Lesser the molecular diffusivity larger the actual concentration in the flow region. An increase in the chemical reaction 
parameter 'K' results in a depreciation in the concentration in the entire flow region. 

5. Higher the Lorentz force larger the Nusselt number at both the walls. An increase in the radiation absorption parameter 
Qi enhances INul at y=+l and reduces at y = -1. 

6. The rate of mass transfer enhances with Schmidt number Sc. The rate of mass transfer enhances with K at y = +1 while 
at y = -1 it enhances with K < 1.5 and reduces with higher K > 2.5. 
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